Wavelet analysis of epileptic spikes 
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Interictal spikes and sharp waves in human EEG are characteristic signatures of epilepsy. These 
potentials originate as a result of synchronous, pathological discharge of many neurons. The reliable 
detection of such potentials has been the long standing problem in EEG analysis, especially after 
long-term monitoring became common in investigation of epileptic patients. The traditional defini- 
tion of a spike is based on its amplitude, duration, sharpness, and emergence from its background. 
However, spike detection systems built solely around this definition are not reliable due to the pres- 
ence of numerous transients and artifacts. We use wavelet transform to analyze the properties of 
EEG manifestations of epilepsy. We demonstrate that the behavior of wavelet transform of epilep- 
tic spikes across scales can constitute the foundation of a relatively simple yet effective detection 
algorithm. 

PACS numbers: 87.10. +e, 05.45.-a 



Recordings of human brain electrical activity (EEG) 
have been the fundamental tool for studying the dynam- 
ics of cortical neurons since 1929. Even though the gist 
of this technique has essentially remained the same, the 
methods of EEG data analysis have profoundly evolved 
during the last two decades. In 1985 Babloyantz et al. 
demonstrated that certain nonlinear measures, first in- 
troduced in the context of chaotic dynamical systems, 
changed during slow- wave sleep Q. The flurry of re- 
search work that followed this discovery focused on the 
application of nonlinear dynamics in quantifying brain 
electrical activity during different mental states, sleep 
stages, and under the influence of the epileptic process 
(for a review see for example 0, 01) ■ It must be em- 
phasized that a straightforward interpretation of neu- 
ral dynamics in terms of such nonlinear measures as the 
largest Lyapunov exponent or the correlation dimension 
is not possible since most biological time series, such as 
EEG, are nonstationary and consequently do not sat- 
isfy the assumptions of the underlying theory. On the 
other hand, traditional power spectral methods are also 
based on quite restrictive assumptions but nevertheless 
have turned out to be successful in some areas of EEG 
analysis. Despite these technical difficulties, the number 
of applications of nonlinear time series analysis has been 
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growing steadily and now includes the characterization of 
encephalopaties monitoring of anesthesia depth 
characteristics of seizure activity @, and prediction of 
epileptic seizures Q. Several other approaches are also 
used to elucidate the nature of electrical activity of the 
human brain ranging from coherence measures 0, 3 and 
methods of nonequilibrium statistical mechanics [10| to 
complexity measures [ill Il2| . 

One of the most important challenges of EEG analy- 
sis has been the quantification of the manifestations of 
epilepsy. The main goal is to establish a correlation be- 
tween the EEG and clinical or pharmacological condi- 
tions. One of the possible approaches is based on the 
properties of the interictal EEG (electrical activity mea- 
sured between seizures) which typically consists of lin- 
ear stochastic background fluctuations interspersed with 
transient nonlinear spikes or sharp waves. These tran- 
sient potentials originate as a result of a simultaneous 
pathological discharge of neurons within a volume of at 
least several mm 3 . 

The traditional definition of a spike is based on its 
amplitude, duration, sharpness and emergence from its 
background 0,0]. However, automatic epileptic spike 
detection systems based on this direct approach suffer 
from false detections in the presence of numerous types 
of artifacts and non-epileptic transients. This shortcom- 
ing is particularly acute for long-term EEG monitoring 
of epileptic patients which became common in 1980s. To 
reduce false detections Gotman and Wang 15] made the 
process of spike identification dependent upon the state 
of EEG (active wakefulness, quiet wakefulness, desyn- 
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chronized EEG, phasic EEG, and slow EEG). This modi- 
fication leads to significant overall improvement provided 
that state classification is correct. 

Diambra and Malta ^(| adopted nonlinear prediction 
for epileptic spike detection. They demonstrated that 
when the model's parameters are adjusted during the 
"learning" phase to assure good predictive performance 
for stochastic background fluctuations, the appearance 
of an interictal spike is marked by a very large forecast- 
ing error. This novel approach is appealing because it 
makes use of changes in EEG dynamics. One expects 
good nonlinear predictive performance when the dynam- 
ics of the EEG interval used for building up the model is 
similar to the dynamics of the interval used for testing. 
However, it is uncertain at this point whether it is possi- 
ble to develop a robust spike detection algorithm based 
solely on this idea. 

As Clark et al. put it succinctly, automatic EEG anal- 
ysis is a formidable task because of the lack of "...features 
that reflect the relevant information "^3- Another dif- 
ficulty is the nonstationary nature of the spikes and the 
background in which they are embedded. One technique 
developed for the treatment of such nonstationary time 
series is wavelet analysis 0, ^| . The goal of this paper 
is to characterize the epileptic spikes and sharp waves in 
terms of the properties of their wavelet transforms. In 
particular, we search for features which could be impor- 
tant in the detection of epileptic events. 

The wavelet transform is an integral transform for 
which the set of basis functions, known as wavelets, are 
well localized both in time and frequency. Moreover, the 
wavelet basis can be constructed from a single function 
ip(t) by means of translation and dilation: 



4>a;; 



1> 



t-t 



(1) 



ip(t) is commonly referred to as the mother function or 
analyzing wavelet. The wavelet transform of function 
h(t) is defined as 
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W(a,t ) = ^ \ 
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where denotes the complex conjugate of ip(t). The 

continuous wavelet transform of a discrete time series 
{fti}^ 1 of length N and equal spacing St is defined as 



5t N ^ 



W n (a) = V — > h n 'i> 



n'=0 



(ri - n)St 



(3) 



The above convolution can be evaluated for any of N 
values of the time index n. However, by choosing all N 
successive time index values, the convolution theorem al- 
lows us to calculate all N convolutions simultaneously in 
Fourier space using a discrete Fourier transform (DFT). 
The DFT of {MiLo* is 
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where k — 0, N — 1 is the frequency index. If one 
notes that the Fourier transform of a function ip(t/a) is 
\a\\j}{af) then by the convolution theorem 




W n (a) = VaStyh n ij*(af k ) 



(5) 



frequencies fk are defined in the conventional way. Us- 
ing JSJ) and a standard fast Fourier transform (FFT) rou- 
tine it is possible to efficiently calculate the continuous 
wavelet transform (for a given scale a) at all n simul- 
taneously [2(J. It should be emphasized that formally 
equation JSJ does not yield the discrete linear convolu- 
tion corresponding to @ but rather a discrete circular 
convolution in which the shift n' — n is taken modulo 
N. However, in the context of this work, this problem 
does not give rise to any numerical difficulties. This is 
because, for purely practical reasons, the beginning and 
the end of the analyzed part of data stream are not taken 
into account during the EEG spike detection. 

From a plethora of available mother wavelets, we em- 
ploy the Mexican hat 
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which is particularly suitable for studying epileptic 
events. 

In the top panel of Fig. ^ we present two pieces of 
the EEG recording joined at approximately t = Is. The 
digital 19 channel recording sampled at 240 Hz was ob- 
tained from a juvenile epileptic patient according to the 
international 10-20 standard with the reference average 
electrode. The epileptic spike in this figure (marked by 
the arrow) is followed by two artifacts. The bottom panel 
of Fig. n displays the contour map of the absolute value 
of Mexican hat wavelet coefficients W(a, to). It is ap- 
parent that the red prominent ridges correspond to the 
position of either spike or the motion artifacts. What 
is most important, for small scales, a, the values of the 
wavelet coefficients for the spike's ridge are much larger 
than those for the artifacts. The peak value along the 
spike ridge corresponds to a = 7. In sharp contrast, for 
the range of scales used in Fig. ^ the absolute value of 
coefficients W(a, to) for the artifacts grow monotonically 
with a. 

The question arises as to whether the behavior of the 
wavelet transform as a function of scale can be used to 
develop a reliable detection algorithm. The first step in 
this direction is to use the normalized wavelet power 



w(a,t ) = W 2 (a,t )/a 2 
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instead of the wavelet coefficients to reduce the depen- 
dence on the amplitude of the EEG recording. In the 
above formula a 2 is the variance of the portion of the sig- 
nal being analyzed (typically we use pieces of length 1024 
for EEG tracings sampled at 240 Hz). In actual numeri- 
cal calculations we prefer to use the square of w(a, to) to 
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merely increase the range of values analyzed during the 
spike detection process. In Fig. [21 u> 2 for the signal used 
in Fig. \I\is plotted for three scales A — 3, B — 7 and 
C = 20. 

In the most straightforward approach, we identify an 
EEG transient potential as a simple or isolated epileptic 
spike if and only if: 

• the value of w at a = 7 is greater than a predeter- 
mined threshold value T\, 

• the square of normalized wavelet power decreases 
from scale a — 7 to a = 20, 

• the value of w 2 at a = 3 is greater than a predeter- 
mined threshold value 

The threshold values T\ and T2 may be considered as 
the model's parameters which can be adjusted to achieve 
the desired sensitivity (the ratio of detected epileptic 
events to the total number of epileptic events present 
in the analyzed EEG tracing) and selectivity (the ratio 
of epileptic events to the total number of events marked 
by the algorithm as epileptic spikes). 

While this simple algorithm is quite effective for simple 
spikes such as one shown in Fig. Q]it fails for the com- 
mon case of an epileptic spike accompanied by a slow 
wave with comparable amplitude. The example of such 
complex is given in Fig. E{a). The overlap of the neg- 
ative tail of the Mexican hat with the slow wave yields 
the inherently low values of w 2 at scale .4 (panel (1))) and 
scale B (panel (c)) as compared to those characteristic of 
the "isolated" spike. Nevertheless, the normalized wavelet 
power does decrease from scale B to C. Consequently, 
in the same vein as the argument we presented above, 
we can develop an algorithm which detects the epileptic 
spike in the vicinity of a slow wave by calculating the 
following linear combination of wavelet transforms: 

W{a,t ) = c 1 W(a,t )+c 2 W(a s ,t +T) (8) 

and checking whether the square of corresponding nor- 
malized power w(a,t ) = W 2 (a, t )/a 2 at scales a = 3 
and a — 7 exceeds the threshold value T\ and T 2 , respec- 
tively. The second term in JHJ allows us to detect the 
slow wave which follows the spike. The parameters a s 
and t are chosen to maximize the overlap of the wavelet 
with the slow wave. For the Mexican hat we use a s = 28 
and t = 0.125s. By varying the values of coefficients ci 
and C2, it is possible to control the relative contribution 
of the spike and the slow wave to the linear combination 



For testing purposes, we built up the database of arti- 
facts and spikes. We made available some of these EEG 
tracings at [2l| along with the examples of the numeri- 
cal calculations. While the analysis of the pieces of EEG 
recordings such as those shown in Fig. [3]and[3]is essential 
in determining the generic properties of epileptic events, 
it can hardly reflect the difficulties one can encounter 
in interpretation of clinical EEG. Therefore we selected 
four challenging EEG tracings with 340 epileptic events. 
The algorithm described in this work marked 356 events 
out of which 239 turned out to be the epileptic events. 
Thus the sensitivity of the algorithm was 70% and its se- 
lectivity was equal to 67%. We then analyzed the same 
tracings with the leading commercial spike detector de- 
veloped by the Persyst Development Corporation (In- 
sight 2001.07.12). This software marked 654 events out 
of which 268 were epileptic events. Thus slightly better 
sensitivity of 79% was achieved at the expense of the low 
41% selectivity. The performance of preliminary numeri- 
cal implementation of the detection algorithm presented 
in this work is excellent and allows to process 24 hour 
EEG recording (19 channels sampled at 240 Hz) in a 
matter of minutes on the average personal computer. 

The goal of wavelet analysis of the two types of spikes, 
presented in this paper, was to elucidate the approach 
to epileptic events detection which explicitly hinges on 
the behavior of wavelet power spectrum of EEG signal 
across scales and not merely on its values. Thus, this 
approach is distinct not only from the detection algo- 
rithms based upon discrete multiresolution representa- 
tions of EEG recordings [H |H |H El 13 but also from 
the method developed by Senhadji and Wendling which 
employs continuous wavelet transform p7) |. 

Epilepsy is a common disease which affects 1-2% of 
the population and about 4% of children |2§|- In some 
epilepsy syndromes interictal paroxysmal discharges of 
cerebral neurons reflect the severity of the epileptic dis- 
order and themselves are believed to contribute to the 
progressive disturbances in cerebral functions ( eg. speech 
impairment, behavioral disturbances) |29j . In such cases 
precise quantitative spike analysis would be extremely 
important. The epileptic event detector described in this 
paper was developed with this particular goal in mind 
and its application to the studies of the Landau-Kleffncr 
syndrome will be presented elsewhere. 
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FIG. 1: Top panel: simple epileptic spike (marked by S) fol- 
lowed by two artifacts. Bottom panel: contour map of the 
absolute value of the Mexican hat wavelet coefficients (arbi- 
trary units) calculated for the EEG signal shown above. The 
shades of blue correspond to low values and the shades of red 
to high values. 




FIG. 2: Square of normalized wavelet power for three different 
scales A < B < C (Panels (b)-(d)). The EEG signal shown 
in panel (a) is the same as the one used in Fig. 
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FIG. 3: (a) Epileptic spike - slow wave complex. The ampli- 
tude of the slow wave is comparable to that of the spike. The 
square of normalized wavelet power for this signal is shown in 
panels (b)-(d) for three different scales A < B < C. 



